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'■ ABSTRACT 

C*) ^, Context. Cosmic rays are confined in the atmospheres of galaxy clusters and, therefore, they can play a crucial role in the heating 
I ■ of their cool cores. 

O ' Aims. We discuss here the thermal and non-thermal features of a model of cosmic ray heating of clusters' cores that can provide a 
, solution to the cooling-flow problems. To this aim, we generalize a model originally proposed by Colafrancesco, Dar & DeRujula 
^ ' (2004) and we show that our model predicts specific correlations between the thermal and non-thermal properties of galaxy 
clusters and enables various observational tests. 

Methods. The model reproduces the observed temperature distribution in clusters by using an energy balance condition in which 
the X-ray energy emitted by clusters is supplied, in a quasi-steady state, by the hadronic cosmic rays, which act as "warming 
^ , rays" (WRs). The temperature profile of the intracluster (IC) gas is strictly correlated with the pressure distribution of the 
WRs and, consequently, with the non-thermal emission (radio, hard X-ray and gamma-ray) induced by the interaction of the 
J:^ ■ ■ WRs with the IC gas and the IC magnetic field. 

Results. The temperature distribution of the IC gas in both cool-core and non cool-core clusters is successfully predicted from 
the measured IC plasma density distribution. Under this contraint, the WR model is also able to reproduce the thermal and 
non-thermal pressure distribution in clusters, as well as their radial entropy distribution, as shown by the analysis of three 
clusters studied in details: Perseus, A2199 and Hydra. The WR model provides other observable features of galaxy clusters: a 
correlation of the pressure ratio (WRs to thermal IC gas) with the inner cluster temperature {Pwr/ Pth) ~ {kTi„„er)~^^^ , a 
correlation of the gamma-ray luminosity with the inner cluster temperature L~f ~ {kTi„„er)*^^ , a substantial number of cool-core 
clusters observable with the GLAST-LAT experiment, a surface brightness of radio halos in cool-core clusters that recovers the 
' observed one, a hard X-ray ICS emission from cool-core clusters that is systematically lower than the observed limits and yet 
rN , observable with the next generation high-sensitivity and spatial resolution HXR experiments like Simbol-X. 

• Conclusions. The specific theoretical properties and the multi-frequency distribution of the e.m. signals predicted in the WR 
■ ■ model render it quite different from the other models so far proposed for the heating of clusters' cool-cores. Such differences 

make it possible to prove or disprove our model as an explanation of the cooling-fiow problems on the basis of multi-frequency 
observations of galaxy clusters. 
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1. Introduction 

Theoretical description of the intra-cluster (IC) plasma 
in the central regions of many clusters predicts that it is 
cooler than in their outskirts, because it radiates X-rays at 
such a rate that the plasma cooling time is much shorter 
than the age of the cluster. However, such rapid cooling 
is not observed and the central-temperature depression in 
cluster's cores is not as deep as expected on the basis of the 
plasma's cooling rate: the observed central temperature 
Tinner settlcs, in fact, at a fraction ^ 1/3 — 1/2 of the outer 
temperature Toutcr (see, e.g. McNamara 1997, Peterson 



et al. 2003, Piffaretti et al. 2005, Donahue et al. 2005, 
see also Bregman 2004 for a review). The lack of cool- 
ing gas with T < Tinner in cluster cores is also supported 
by the absence of line emission corresponding to gas be- 
low Tinner, f-s showii by high spatial resolution imaging 
with Chandra (e.g. McNamara et al. 2000; Fabian 2000; 
Blanton et al. 2001; Allen et al. 2001; Lewis ct al. 2002; 
Blanton et al. 2003) and high spectral resolution measure- 
ments with the XMM-Newton RGS instrument (Peterson 
et al. 2001, Tamura et al. 2001, Kaastra et al. 2001; Kahn 
et al. 2002; Peterson et al. 2002) and EPIC instrument 
(Bohringer et al. 2001, 2002; Molendi & Pizzolato 2001; 
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Matsushita et al. 2002). These are referred here as the 
cooUng-flow (CF) problems. 

It has been widely recognized that some form of heat- 
ing is necessary to quench the cooling flow and form a 
warm-core in the IC gas structure. Available sources of 
heating in clusters cores are provided by AGN jets and 
lobes, pressure waves, buoyant bubbles and cavities, intra- 
cluster shock waves, leptonic and hadronic cosmic-rays. 
The presence of a central radio source with jet and lobes, 
the mechanical heating provided by pressure waves pro- 
duced by the jetted AGNs, the formation of buoyant bub- 
bles filled with energetic plasma, the thermal conduc- 
tion from the hot outer layers of the clusters, may cer- 
tainly alleviate the CF problems (see e.g., Bohringer et al. 
2002; Churazov et al. 2002; Ruszkowski & Begelman 2002, 
Fabian 2004, Voit & Donahue 2005, Briiggen & Kaiser 
2002, Ruszkowski et al. 2004, Vernaleo & Reynolds 2006, 
Reynolds et al. 2005). 

Some difficulties with solutions along these lines, however, 
might stand out: a first difficulty is that not all cool- 
ing flow clusters contain a powerful radio source at their 
center (e.g., the cluster RXJ0820. 9+0752; Bayer-Kim et 
al. 2002, and the clusters A1650 and A2244, Donahue et al. 

2005) , nor a large amount of spherically distributed buoy- 
ant bubbles out to the cluster's core boundary. Another 
difficulty is that the bremsstrahlung cooling rate is pro- 
portional to ri^T^/^, with ric the plasma's electron num- 
ber densitjQ and T its local temperature: a suitable con- 
duction and/or heating and pressure-building mechanism 
must somehow adapt itself to this behaviour, particularly 
as a function of no(r), which varies by orders of magnitude 
along the clusters' radii. Other difficulties stand out with 
the spatial distribution of the energy deposited by AGN 
jets and cavities: on one hand, the energy released from 
AGN jets and cavities seems to be sufficient to quench 
cooling flow and increase the central gas entropy (e.g., 
Sijacki & Springel 2006); on the other hand, it is not yet 
clear if such energy can be properly distributed to repro- 
duce the observed cool-core temperature profiles (see, e.g., 
discussion by Vernaleo & Reynolds 2005, 2006, Heinz et 
al. 2006, Vernaleo & Reynolds 2007). SmaU duty cicles 
and a higher fraction of AGNs in cluster cores (see Bird 
et al. 2007) might however weaken the difficulties related 
to the AGN heating scenarios. In addition, there seems to 
be a small AGN activity in groups (Dwarakanath & Nath 

2006) . 

Sound waves and shocks do not seem to be efficient in 
quenching cooling fiows (see, e.g., Fujita & Suzuki 2006). 

Irrespective of the nature of the heating source, it 
seems that a general property of the heating agent - which 
derives from the smoothness and from the similarity of the 
heating distribution required to quench cooling in several 
clusters - is that they must be spatially distributed in the 
whole cluster core. In fact, it has been shown (Fabian & 
Sanders 2006, Fabian et al. 2006, Sanders & Fabian 2007) 

^ Hydrogen in the plasma is fully ionized, and we are using 
here rip « ne. 



that the non-thermal pressure in the Perseus cluster (one 
of the best studied clusters with cool cores and mini radio- 
halos) is quite similar (slightly steeper) to the thermal one. 
In this context, it has been shown that a CR population 
in clusters that is radially distributed as the radial pro- 
file of the IC gas, is able to recover the observed radial 
structure of the IC gas temperature in several clusters 
(Colafrancesco, Dar & DeRujula 2004). The CRs in this 
last theory are not confined to the disk of the Galaxy, but 
permeate a much larger halo, being produced by "cannon- 
balls" travelling for kiloparsecs in the interstellar medium 
of galaxies, decelerating by knocking out its constituents, 
which are thereby accelerated to become CrU (see, e.g., 
DeRujula 2004 for a review, and references therein). 
Beyond the debate on the cannonball picture (see, e.g., 
Hillas 2006; Dar & DeRujula 2006) and its relevance for 
the CR origin, we stress here that any mechanism which 
produces radial distributions of WRs similar to the profile 
of the IC gas can provide quenching of cooling flows by 
the action of WRs heating. 

In this paper we consider and justify, on the basis of a 
consistent picture of CR diffusion in cluster atmospheres 
(see Appendix A), that the WR number density nwnir) 
can be, in fact, distributed with a radial proflle that is 
similar to that of the IC gas, and follows the general form 
nwR ~ n'^, i-e. it is proportional to a power a of the IC 
gas density Ue {r) . The total amount of WRs in our model 
and the slope a of the nwRir) profile are left as free pa- 
rameters that will be subsequently constrained by fitting 
the temperature profile of the IC gas in each specific clus- 
ter core. The resulting CR density profile will be further 
tested against the entropy and pressure profiles in the clus- 
ter cores, the spectrum and surface brightness profile of 
the cluster radio halos and the future observations in the 
hard X-ray and gamma-ray bands. 

We will first show that the presence of WRs in cluster cores 
are able, in fact, to heat the IC gas and reproduce their 
radial temperature profiles. We will then start from this 
constraint to discuss the consequences that such a popu- 
lation of WR hadrons (protons), produced and stored in 
the atmospheres of galaxy clusters, has for the gamma- 
ray emission, for the hard X-ray (HXR) emission and for 
the presence of diffuse radio emission of both cool-core 
and non cool-core clusters. We will finally show that it 
is possible to use these predictions, in the framework of a 
multi- wavelength observational strategy, to set constraints 
on the amount and spatial distribution of WRs in clus- 
ters, and therefore disentangle between the WR model 
discussed here and other models for the onset of warm- 
cores in galaxy clusters. 

The structure of the paper is the following: in Sect. 2 we 
delineate the theoretical framework that is able to model 
the temperature structure of the cluster cool cores and 
determine the main parameters (i.e., the radial distribu- 
tion and the density of WRs) that fit the cluster tempera- 

^ When referring to CRs in CF clusters, we use the expres- 
sions "CRs" and "WRs" as synonyms. 
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ture profiles. In Sect. 3 we will also discuss the density and 
pressure distribution of WRs in the clusters' atmospheres 
that is consistent with the temperature profile. In Sect. 4 
we will apply our procedure to a sample of ten well stud- 
ied clusters. These clusters have been chosen according to 
the following criteria: i) 4 clusters which have cool cores 
and diffuse radio emission (these clusters certainly have 
a population of CRs in their centers); ii) 2 clusters with 
radio halos (i.e. with CRs in their atmospheres) but with- 
out cool cores, in order to compare their properties with 
the radio-active, cool-core clusters; iii) 4 clusters with cool 
cores but with no evident central diffuse radio emission 
in order to compare these non radio-active clusters with 
the previous ones. In Sect. 4 we also study the thermal and 
non-thermal pressure structure of the selected clusters and 
we will show that the pressure ratio Pwr/ Pth found at the 
cluster center correlates with the central gas temperature 
Tinner providing evidence for the basic regulation mecha- 
nism induced by WRs. We will discuss the predictions of 
our model for some specific clusters (like Perseus, A2199 
and Hydra) for which detailed thermal and non-thermal 
pressure information is available as well as information 
on their radial entropy distribution. In Sect. 5 we will dis- 
cuss the expected gamma-ray luminosity and flux of the 
cluster sample here considered that is correlated with the 
cool-core temperature structure predicted in our model. In 
Sects. 6 and 7 we will present analogous considerations for 
the expected radio and HXR emission that are constrained 
by the cluster temperature structure. We will discuss the 
differences between our model and other scenarios for the 
cool-core heating in Sect. 8 and we will resume our main 
conclusions in Sect. 9. 

Throughout the paper, we use a flat, vacuum-dominated 
cosmological model with — 0.3, JIa = 0.7 and h = 0.7. 

2. Modelling cluster cool cores 

The virial expectation for the temperature of an un- 
magnetized plasma in a spherical cluster with mass M(r) 
within a radius r, is: 



fcT(r) 



3r 



2keV 



M(r) 



1013 M, 







100 kpc 



(1) 



(see, e.g. Colafrancesco, Dar & De Rujula 2004). The ob- 
served outer temperature Toutcr of the IC plasma - at large 
radii where the cooling time should be longer than the 
cluster's age - is roughly compatible with the virial expec- 
tation. In the inner regions of many clusters the tempera- 
ture is somewhat smaller (~ 1/2 — 1/3 of Touter) than that 
obtained form Eq. ([1]), but signiflcantly larger than the 
cooling rate would imply (see, e.g., Peterson et al. 2001). 
To reproduce the observed temperature proflle from the 
inner cool region to the outer hot region, we will consider 
a population of WRs spatially distributed in the cluster 
atmosphere (as obtained by solving a transport equation, 
see Appendix A) that efficiently heats the IC gas. 



The effect of the "warming rays" on a cluster's tem- 
perature distribution will be computed under the following 
simplifying assumptions: 

i) spherical symmetry; 

ii) the radial distribution of the WR density is approxi- 
mately proportional to a power a of the thermal plasma 
density, i.e. nwR oc n" , with a and the normalization con- 
stant being free parameters. In the Appendix A we show 
that the effect of a simple model of cosmic ray diffusion in 
cluster atmospheres can provide such kind of radial distri- 
bution, for quite general CR source terms; 

iii) the temperature of the plasma in its "initial state" 
(as the cluster is born) is approximated by a radius- 
independent constant value; 

iv) cooling by X-ray emission and heating and pressure- 
building by WRs are the dominant evolutionary agents. 
Flows and cluster ageing are relatively unimportant. 

Given this framework, the rate per unit volume 
at which the IC plasma loses energy by thermal 
bremsstrahlung emission of X-rays writes as: 

de 
dt 



33 h ml c3 



Gi 



4.8 X lO^^'if- 



1 erg cm^ 



(2) 



where z is an average charge of the IC plasma (we have 
approximated here the Gaunt factor G by unity) . We write 
the IC gas number density as ne(r) = rieo ■ gif), where the 
function g{r) contains all the radial dependence of the IC 
gas density and is obtained from X-ray observations of 
each specific cluster. The rate per unit volume at which 
the WRs deposit energy in the IC plasma is proportional 
to nwRir) ■ ne{r) and it is therefore: 



de 



= ^ •["WR('')]-K(f)] = fe-"WR,0-"cO-[5'(f)]^ 



(3) 



since we have assumed that the WR number density writes 
as nwR = nwR,o ■ [9(1^)]°'^ i-e. proportional to [no(r)]". 
The evolution of the temperature anywhere in the IC 
plasma satisfies the energy-conservation relation: 



3 k nc{r) 



dT{r, t) 
dt 



de 
'di 



de 
dt 



(4) 



with the initial condition r(r, 0) — Ti. We can rewrite 
Eq. (gl) as: 

dT{r,t) a?io(r) 



VT- {b/a)[nwR{r)/ne{r)] 
and integrate, to obtain: 



3fc 



■dt 



(5) 



In 



^T{r,t)/T,-1 



lT{r, t) 



V Tf 

bnwRir)nc{r) 
QnJr)kTf 




{t - ti) 



(6) 
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where 



kT[ = [{b/a){nwR,o/neo)g{r)°'' 



(7) 



is the steady temperature first reached in the central re- 
gion of the cluster at small radii r , where the density ne{r) 
is high, and at asymptotically late time in the periphery, 
where the density ne(r) is vanishingly small. The r.h.s. 
of Eq. © has been written in a way that shows how the 
evolution of T{r, t) is determined by the ratio between the 
WR deposition rate, 6nwR(?')ne(r), and the thermal en- 
ergy per unit volume, ne{r)kTf, in the inner region where 
the cluster's temperature is Tf. 

The temperature evolution for a cluster similar to 
A2199 (we assume here indeed the same IC gas density 
profile of A2199 and its outer temperature) as described 
by Eq. is shown in Fig. [T] at epoch t = 0, t^/3, and 
3*,^, where is the current age of the Universe. This fig- 
ure shows that, in our model, a cluster starts as an isother- 
mal system at high redshifts (i.e., t ^ 0) and then a warm- 
core develops with time under the influence of cooling and 
heating effects balancing each other in order to provide a 
stationary configuration at each epoch. We note that the 
specific choice of the initial epoch ti is not very important 
in the solution of eq.®, because the epoch t (redshift z) at 
which the cluster is observed is much larger (smaller) than 
the initial epoch ti (redshift Zi ^ z). We have nonetheless 
calculated the temperature evolution by assuming differ- 
ent values of the initial redshift, 50 < z{ti) < 200, and we 
found that no appreciable difference appears in the final 
state of the temperature radial distribution. 

In cool-core clusters the central density is high enough 
for the equilibrium between X-ray cooling and WR heating 
to have been reached at the cluster's center in a time of 
the order of its current age. In such a case, Tf ~ Tinner, i-e. 
the presently observed inner-cluster temperature. For all 
the clusters we study, the cooling time in the outer parts 
of the cluster is much longer than t^, so that T; ~ Touter- 
According to the value of Ti, of the IC gas density Ue 
and of the time t, Eq.® predicts that some clusters have 
reached this stage, while others have not. 

The procedure here described is used to fit the tem- 
perature profile of clusters with and without cool cores. 
Given the bremsstrahlung energy loss rate, the tempera- 
ture radial profile of each cluster depends on the density 
and on the spectrum nwFt.{E,r) of WRs in clusters. 

3. Warming Rays in cool core clusters 

We consider a spectrum of relativistic WRs (protons) as 
given by: 



NwB{E,r) = NwRflE-'[g(r)]'' , 



(8) 



with s typically in the range ~ 2.3 — 3.3 (see, e.g., 
Marchegiani et al. 2007), and we use s — 2.7 as a ref- 
erence value, hereafter, for energies larger than Emin = 
rupc^ ■ [1 + 3.4 • lQ-^{kT/keV)] (see, e.g. Furlanetto & 
Loeb 2002). The value of E^ax is usually very large and 
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Fig. 1. The evolution of the IC gas temperature profile pre- 
dicted for a cluster similar to A2199 is shown at times t = Q 
(solid), tjj /i (dashed), t^j (dotted) and Si^^ (dot-dashed), with 
tjj being the age of the Universe. We assume an initial temper- 
ature of 4.2 keV and a density profile equal to that observed 
in A2199. Values a — 0.83 and nwR,o = 2.3 • 10"'^ cm~^ are 
assumed for this specific calculation. 

is therefore irrelevant in our calculations for spectral in- 
dices s > 2. 

The quantity N\yR,o can be derived by assuming that 
the WR heating rate b ■ nwR • provided by the energy 
lost by the WRs when they interact with the IC medium, 
recovers the IC gas temperature profile through eq.®. 

The fundamental processes by which the warming 
rays can heat the IC gas are mainly Coulomb and 
hadronic interactions with the plasma nuclei and mag- 
neto hydro-dynamical (MHD) waves excited by CR proton 
anisotropics. 

For moderately relativistic CR nuclei with a charge z, 
velocity v = P c and Lorentz factor 7, the Coulomb energy- 
loss rate is: 



fa K z 

dt 



In 



2 mo P^j"^ 



(9) 



where is the (suitably averaged) squared charge of the 



plasma's nuclei, K 
2.82 fm. Here /„ 



4 Trno rric 
hujp, with Up — 



and 



'/me 



[47rne e^/me]^/^ the 
plasma frequency!! The heating rate given in Eq. ([9]) , for 
typical plasma densities near the center of clusters {uc ^ 



It has been noticed by Colafrancesco, Dar & DeRujula 
(2004) that a highly ionized plasma is much more efficient than 
non-ionized gas in slowing down the CR nuclei. Energy losses 
via Coulomb collisions with neutral atoms become significant 
only for energy transfers larger than their ionization potential 
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10~^ cm~'^), is ~4 times larger than for neutral hydrogen. 
For protons (2 = 1) with a typical WR energy (7 = 2) and 
for a plasma consisting of ~ 0.93% H, ~ 0.07% He (by 
number) and traces of heavier elements, one finds Z"^ ~ 
1.14 and dEjdt k, 10"^'' GeV/s, for the quoted electron 
number density. 

Hadronic collisions between WRs and plasma nuclei 
constitute an energy-loss mechanism which is approxi- 
mately as efficient as that given in Eq. ([9]). Here we con- 
sider the collisions of WR protons with ambient protons, 
the discussion of proton-nucleon and nucleon-nucleon col- 
lisions being far too cumbersome to be justified by the 
small admixtures of nuclei heavier than protons in the 
WR and IC constituencies. For values 7 > 2 of the inci- 
dent proton, the cross section is Cpp ~ 4 • 10~^^ cm^, 
and it is dominantly inelastic. In practically all of these 
collisions, the incident proton survives unscathed, as the 
"leading" (most energetic) final-state particle, carrying on 
average ~ 70% of the incident proton's energy. Thus, a 
proton's energy-loss rate by hadronic collisions in a hy- 
drogenic plasma is 



dE 

—— ~ 0.7 n„ cr„„ c E 
dt 



(10) 



(Mannheim & Schhckeiser 1994), that reads dE/dt^l.6 x 
10^^^ GeV/s, for our reference value of and 7 = 2. 
This is very close to the Coulomb energy loss previously 
computed for the same reference parameters. 

An additional process which has been discussed in the 
literature is the effect of the MHD waves excited by proton 
anisotropy, which might occur as relativistic, charged par- 
ticle stream faster than the local Alfven speed (see, e.g. 
Skilhng 1971); these waves can be dissipated efficiently 
in the background ICM via Landau damping and other 
dissipative mechanisms (see, e.g., Foote & Kulsrud 1979). 
The effect of this mechanism in cluster atmospheres has 
been studied in details by Loewenstein et al. (1991): under 
the assumptions that the wave velocity is large enough to 
heat the IC gas significantly, and that the Coulomb and 
hadronic heating can be neglected, these authors found 
that the thermal pressure becomes overwhelmed by the 
CR proton pressure. Therefore, they concluded that such 
CR heating mechanism alone is not able to stabilize the 
cooling flows. 

Following the previous considerations, we assume 
therefore that the Coulomb and hadronic heating are the 
dominant sources and we search for the amount and spa- 
tial distribution of WRs required to produce the temper- 
ature structure observed in a sample of clusters with cool 
cores (i.e. where heating is required) and radio halos (i.e., 
where CRs are required). In Sect. 14.31 we will check our 
results by comparing the non-thermal and thermal pres- 
sures predicted in our model with those observed in the 
same clusters. 



I. The corresponding small impact parameters imply relatively 
small cross sections. For the IC plasma, the role of I is played 
by the much smaller quantity /p. 



The total WR loss rate is given by the sum of eqs. 
and (Uni) 



dE\ 

/ tot 



(11) 



with 



Ahadr — 0.7 appCE 



(12) 



and is approximately proportional to He (by neglecting the 
logarithmic dependence on rig contained in ojp in eql9]). 

The value of b (see eqI3]) is found by convolving the 
energy loss rate with the WR spectrum and equating with 
the warming rate 



de 
dt 



NwR(E,r) 



E„ 



dE 
'dt 



dE 



WR 

which writes as 



b ■ nwR{r) ■ ne{r) = NwR,og°'{r)ne{r)lA , 
where 

Ia= E~^{Acoul + Ahadr)dE . 



(13) 



(14) 



(15) 



Since the WR density at the cluster's center is given by 
nwR,o = / dENwRfiE-" , (16) 



E„ 



we obtain 



Ia 



Je 



dEE- 



(17) 



Inserting the previous expression of b in eq.® we obtain 
the value of nwR,o required to fit the temperature profile, 
and in turn the pressure of the WRs as 

1 n^"- 

PwRir) = ^mpCj NwRi'-f,r)f3'yv{'-f)dj , (18) 
where 7 = Ep/{mp(?). 



4. The effect of WRs in cluster cores 

4.1. Temperature profiles 

The solution of eq.® provides a fit to the cluster's tem- 
perature profile by adjusting the free parameters a and 
nwRfi- We have applied this analysis to ten well stud- 
ied clusters, eight of which are cool-core (CC) clusters 
and two are non cool-core (NCC) clusters (see Table [1]). 
Specifically, we considered 4 clusters with cool-cores which 
have central radio emission (usually of the class of mini 
radio-halos) : these are A2199, Perseus, A2390 and RX 
J1347. 5-1145. For comparison, we also considered two well 
known non cool-core clusters with radio halos (Coma and 
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Cluster 


z 


Tinner 


Touter 


Notes 


Ref. 






keV 


keV 






A262 


0.0162 


0.95 


2.4 


CC 


[1] 


A2199 


0.030 


1.6 


4.2 


CC; RH 


[2] 


A133 


0.057 


2.05 


4.61 


CC; Relic 


[1] 


Perseus 


0.0179 


3.04 


6.7 


CC; RH 


[3] 


Hydra 


0.0539 


3.1 


4.0 


CC 


[4] 


A 1795 


0.0625 


3.5 


6.5 


CC 


[1] 


A2390 


0.2304 


4.89 


9.81 


CC; RH 


[1] 


RX J1347.5-1145 


0.451 


7.04 


19.4 


CC; RH 


[5] 


Coma 


0.023 


8.2 


8.2 


NCC; RH 


[6] 


A2163 


0.203 


14.6 


14.6 


NCC; RH 


[7] 



Table 1. The list of clusters analyzed. Redshift z, inner observed temperature Tinner and outer observed temperature Toutcr 
are given. The last two columns show which cluster has a radio halo (CC= cool-core; NCC= non cool-core; RH= radio halo) 
and the references of the data we used for the fit. [1] Vikhlinin et al. (2005); [2] Johnstone et al. (2002); [3] Churazov et al. 
(2003); [4] David et al. (2001); [5] Allen et al. (2002); [6] Briel et al. (1992); [7] Elbaz et al. (1995). 



A2163) and some of the best studied clusters with cool- 
cores but no evident diffuse radio emission, i.e. A262, 
A133 (which has a relic radio emission feature), Hydra 
and A1795. 

For each cluster in Tab[l] we found the values of a and 
'^WR.o which best fit the observed temperature profile, by 
assuming that initially (i.e. at high z when the cluster 
formed) each cluster was isothermal with temperature T, 
equal to the outer observed temperature, and that the 
evolution of the temperature proceeds according to eq. ^ 
and finishes at the redshift z at which each cluster is ob- 
served. 

The WR density nwR(7') found through this procedure is 
the one required to heat the IC gas at the level shown by 
the temperature profile by using eq.®. 

We notice that the WR density in cool-core clusters 
is required to be flatter (i.e. with values a < 1) than the 
IC gas density to provide the heating rate necessary to 
quench the gas coohng and reproduce the cluster's tem- 
perature profile. On the other hand, clusters without cool- 
cores require a « 1 because they did not yet developed gas 
cooling conditions (see discussion at the end of Sect. 2). 

4.2. Pressure and density of WRs 

The ratio of the central WR pressure Pwb. to the IC gas 
thermal pressure Pth is found by using eq. lfTS)) and the 
expression Pth = 2ne(r)kT(r) (by including the contribu- 
tion of both thermal protons and electrons) . The value of 
the best- fit central pressure ratio PwrI Pth found for each 
cluster considered in this study is reported in Table [H 

The results of our fitting procedure to the cluster tem- 
perature profiles indicate that there is a strong correla- 
tion between the inner cluster temperature Tinner and the 
pressure ratio PwrI Pth at the cluster center (see Fig.[2])- 



This correlation can be well represented with a power-law 
function PwrJ Pth oc {kTinnerY with 5 ~ —2/3. 
Such correlation is due to the fact that in clusters with 
cooler inner temperature, a high IC gas density triggers 
fast gas cooling; therefore, a stronger warming action 
(compared to the thermal energy) is required to quench 
the cooling and to adjust the temperature profile to the 
observed one by the energy-conservation relation in eq. ([6]) . 
The central pressure ratio PwrI Pth that is required in or- 
der to equilibrate gas cooling and then settle the IC gas at 
the inner temperature value Tinner is obtained from eq.l]?]) 
and reads 

P /P 1 fa\ Pwr 
PwrI Pth - T 

Z^/kl inner \ J nw Rfi 
kTinner ^ ^ ^ 

The correlation PwRlPth oc {kTinner)~'^^^ is explained 
by the combined temperature dependence of the term 
oc {kTinner)^^^'^ and of the quantity (f)(£wfl) ~ 
{kTinner)"^^^ , which depends on Tinner through the lower 
limit of integration Eniin{Tinner) in the expressions of b 
and (swr) (see Figl3]). Note that the Tinner dependence 
of b and (ewr) arises only from the spectral properties of 
WRs in our model, while the scaling with {kTinner)~^^'^ 
arises from the condition set in eq.([6|. 

The other correlation, which is naturally related to the 
previous one, is found - as expected - between the central 
WR density nwR,o and the value of kTnner (see Fig|4]). 
In fact, from eq.© we obtain that 

f^j. (b^nwR^\\ ^20) 
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Cluster 


a 


nwR,o 
cm~^ 


PWR/Pth 


cm ^ s 


1 


erg s 


-1 


Fhxr 
erg cm~2 s~^ 


A262 


0.83 


2.20 


X 


10" 


-3 


1.23 


3.89 X 


10" 


-y 


1.43 


X 


10*^ 


3.87 


X 


10" 


-14 


A2199 


0.83 


2.31 


X 


10" 


-3 


0.92 


8.43 X 


10" 


-9 


1.08 


X 


10« 


3.06 


X 


10" 


-13 


A133 


0.84 


4.56 


X 


10" 


-4 


0.77 


7.30 X 


10" 


10 


3.53 


X 


10*2 


6.10 


X 


10" 


-15 


Perseus 


0.91 


4.98 


X 


10" 


-4 


0.54 


2.20 X 


10" 


-8 


9.91 


X 


10*2 


1.59 


X 


10" 


-13 


Hydra 


0.97 


6.24 


X 


10" 


-4 


0.57 


3.46 X 


10" 


-9 


1.49 


X 


10*3 


2.57 


X 


10" 


-14 


A1795 


0.96 


5.55 


X 


10" 


-4 


0.50 


3.17 X 


10" 


-9 


1.86 


X 


10*3 


2.41 


X 


10" 


-14 


A2390 


0.94 


2.21 


X 


10" 


-4 


0.41 


1.41 X 


10" 


10 


1.39 


X 


10*3 


6.17 


X 


10" 


-16 


RX J1347.5-1145 


0.89 


9.20 


X 


10" 


-4 


0.33 


5.07 X 


10" 


10 


2.37 


X 


10** 


8.66 


X 


10" 


-16 


Coma 


1.00 


1.07 


X 


10" 


-5 


0.29 


1.65 X 


10" 


-8 


1.23 


X 


10*3 


1.67 


X 


10" 


-13 


A2163 


1.00 


1.42 


X 


10" 


-5 


0.21 


1.31 X 


10" 


-9 


9.71 


X 


10*3 


4.85 


X 


10" 


-15 



Table 2. Col.l: cluster name; Col. 2: a; Col. 3: nwR,o (in units of cm"3)| Col. 4: the pressure ratio PwR/Pth at the cluster 
center; Col. 5: the gamma-ray flux in units of cm"^ s~^; Col. 6: the gamma-ray luminosity in units of erg s~^. (The 
gamma-ray luminosity and flux are integrated in the 0.1-10 GeV band). Col. 7: the HXR flux Fhxr. in units of erg cm"^ s~^ 
integrated in the 10-50 keV band. 
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10 
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Fig. 2. The correlation of the central pressure ratio Pwr/ Pth 
with the inner temperature of the cool core Tinner in keV 
is given for each cluster in Tab.l. The flt PwR/Pth oc 
{kTi„„er)~^^^ is shown as a solid line. Blue dots refer to cool- 
core clusters and green dots to non cool-core clusters. 



for r ^ 0. Since the WRs are normalized to obtain T, 



Tinner: ^e Can wrltc 



nwR,o 



a\/kTi, 



■ rieo 



Ie 
Ia 



(21) 



where Ie = /J''".'"" dEE~^ and I a is defined in ea. (fT5|) . 
Then the ratio Ie/Ia depends on i?mm and, as a conse- 
quence, on kTinner] we found that for our range of tern- 



"1,0: 



CL 
CL 




0/ 



kT fkeV 



Fig. 3. The contribution of the term {j){ewR) (dashed curve) 
to the correlation of the pressure ratio Pwr/ Pth with the inner 
temperature of the cool core Tinner is shown. The solid curve is 
the convolution of the terms {^){ewR) and {kTinner)~^^^ (see 
text for details). 



peratures this dependence is Ie /Ia oc {kTi, 
conclusion is that 



-1.7 



nwR,a 



^ {kTinner^ 



The 



(22) 



which is confirmed by our results (see Fig. [5]). 
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For Perseus we find that a central ratio PwR/Pth ~ 
0.54 is required in order to reproduce the inner tempera- 
ture profile; such a value is very well consistent with that 
found between the non-thermal and the thermal pressures 
in the analysis of Sanders & Fabian (2007). It is also re- 
markable that our model is able to reproduce the steeper 
spatial profile of the non-thermal pressure w.r.t. the ther- 
mal pressure profile found by these authors in their anal- 
ysis of the central regions of Perseus (see FigjS]) . We note 
also that the non-thermal pressure we find in this cluster 
is lower than the thermal one: in this respect, we can con- 
clude that our model has not the problems shown by sce- 
nario of CR heating provided by waves excited by protons 
anisotropy (see Loewenstein et al. 1991). In our frame- 
work, the non-thermal pressure component found in the 
core of Perseus can be entirely attributed to WRs. These 
results bring our WR heating model to a strong consis- 
tency level with the results of the available X-ray obser- 
vations (e.g. Sanders & Fabian 2007). As a consequence 



Fig. 4. The correlation of the central WR density nwR,o with 
the inner temperature of the cool core Tinner- Blue dots refer 
to cool-core clusters and green dots to non cool-core clusters. 
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Fig. 5. The correlation of the ratio nwR,o/^cO with the inner 
temperature of the cool core Tinner - The dashed line shows the 
trend nwRfi/neo (x (kTi nrter ] 



Fig. 6. The radial profile of the WR pressure (dashed curve) 
of Perseus is compared with that of the thermal gas (solid 
curve). The radial profile of the WR pressure resembles the 
radial profile of the non-thermal pressure found by Sanders & 
Fabian (2007). 



4.3- Pressure profiles: the cases of Perseus and A2199 

We compare the pressure profiles of the (non-thermal) 
WRs and of the thermal components obtained in our 
model by fitting the IC gas temperature profile to the 
data available for Perseus and A2199. 



of the presence of WRs in the core of Perseus cluster, we 
expect a non-thermal HXR emission due to ICS of the 
secondary electrons from the central region of Perseus: we 
will discuss in Sect. 7 the predictions of our model for the 
HXR emission in the 2-10 keV band where the instruments 
on board Chandra and XMM-Netwon are sensitive. 
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We have also studied the energy density profile of the 
WRs and of the thermal IC gas in A2199. The WR energy 
density has been evaluated as 



(23) 



and again depends on the WR spatial distribution that fits 
the cluster temperature profile. The thermal IC gas energy 
density is evaluated as Eth = 2>ne{r)kT{r) (by including 
the contribution of both thermal protons and electrons). 
The comparison of the two energy densities as evaluated 
in our model is shown in the upper panel of FiglTl Our 
model predicts that the WR energy density at the center 
of A2199 is approximately equal to the IC gas thermal 
energy in order to fit the cluster's temperature profile. 
Similarly to the case of Perseus, we find a radial profile 
of the WR energy density that decays faster with radius 
with respect to the thermal energy density. 
It is possible to recover the energy density ratio of A2199 
as presented in the model of Guo & Oh (2007) by notic- 
ing that the minimum energy cutoff chosen by these last 
authors is ~ 706 MeV, i.e. much higher than the value 
of Emin assumed in this paper (see Sect. 3). In fact, when 
we cut the WR spectrum at this higher value of Emin, 
we recover an energy density ratio which is quite close to 
the value found in Guo & Oh (2007). We verified that 
the residual difference in this ratio is due to the marginal 
effect of the thermal conduction in their model. 

4.4. Entropy profiles 

The WR heating in the cluster core can be also respon- 
sible for the fiattening of the IC gas entropy profile as 
observed in cool-core clusters. We have compared the en- 
tropy profiles predicted in our model for Perseus, A2199 
and Hydra to the available data for these clusters (see 
e.g., Churazov et al. 2003, Cafli et al. 2004, David et al. 
2001, respectively) and we found a very good agreement 
(see FigH]). This is expected, indeed, because our model 
is able to reproduce the cluster temperature profile given 
the observed density profile of the cluster and therefore 
their combination S = kTinner/nV^ which is referred to 
as the IC gas entropy. This analysis bring further robust- 
ness to the ability of our model to explain the structure of 
clusters' cores in terms of a balance between WR heating 
and IC gas cooling. 

5. Gamma-ray emission of cool core clusters 

The distribution of WRs that heat the IC gas to recover 
the observed temperature profiles in their inner regions 
inevitably produces a gamma-ray emission due to vari- 
ous emission mechanisms: i) — *■ tt" — > 7 -I- 7; ii) inverse 
Compton scattering (ICS) due to secondary electrons pro- 
duced by pp TT^ ^ e^] iii) non-thermal bremsstrahlung 
emission of the same secondary electrons. 

From the WR spectra discussed in this paper (see 
Sect. 4), we calculated the gamma-ray emission produced 
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Fig. 7. The radial profile of the WR energy density of A2199 
is compared with that of the thermal gas. The two panels refer 



to i) 

(see, 
Ek 



a proton spectrum with 7min = 1-1- 3.4 x 10 (kT/keV) 
e.g., Furlanetto & Loeb 2002); ii) a proton spectrum with 
= 706 MeV. 



tn,rmn 



by the decay of the tt^ —^ 77 produced in the interac- 
tions between relativistic protons and the thermal gas 
target protons (see, e.g., Blasi & Colafrancesco 1999; 
Marchegiani et al. 2007 for the details of the calculations). 
The contribution of secondary electrons to the gamma- 
ray emission via ICS of CMB photons and non-thermal 
bremsstrahlung is negligible w.r.t. the 7r° 77 fiux (see 
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Fig. 8. The radial profile of the IC gas entropy S of Perseus 
(upper panel), A2199 (mid panel) and Hydra (lower panel) as 
predicted in the WR model. The curves reproduce the observed 
entropy profiles of the three clusters as derived by Churazov et 
al. (2003) for Perseus, by Caffl et al. (2004) for A2199 and by 
David et al. (2001) for Hydra. 



also Marchegiani et al. 2007) and, therefore, we do not 
consider them in the following discussion. 

Note that, at variance with previous works that consid- 
ered the non-thermal (radio, HXR and gamma-ray) emis- 
sion of CRs in clusters' atmospheres (e.g., Colafrancesco 
& Blasi 1998, Blasi & Colafrancesco 1999, Miniati et al. 
2001a, 2001b, Pfrommer & Ensshn 2004, Marchegiani, 
Perola & Colafrancesco 2007, Pfrommer et al. 2007), 
the density of the WRs is not a free parameter in our 
model, because it is fixed by the procedure presented in 



Sect. 2 used to fit the temperature profile of the cluster. 
Therefore, the cluster's gamma-ray emission is univocally 
determined by this condition and it is strongly related to 
the density and temperature structure of the cluster's cool 
core. 

The total gamma-ray fiux and luminosity consistent 
with the cool-core temperature profile of the considered 
clusters are reported in Table [21 The cluster gamma- 
ray luminosity predicted in our model shows an over- 
all increase with the value of Tinner even though with 
some scatter. The total cluster gamma-ray luminosity 



due to pp 



77 decay scales with the den- 



sity of the WRs and of the IC gas according to Lj cx 
nwR{f)n-p{r)R^ , and it is integrated out to the radius at 
which we see target thermal protons from the clusters' X- 
ray images. The condition imposed by the WR heating 
model to fit the inner temperature profile yields a scal- 
ing cx V kTinner {o-/b)n^i:)R^ which can be rewritten as 



^7 CK \J kTinjier 



{a/h)Lx I {kTY^"^ and yields a scaling 



cx 



{kTinner)^ {a /b) by using the observed Lx ~ relation 
(see Arnaud 2005 for a review) and a value kTinner ~ 
(1/3 - l/2)fcr. The scaling of L-y with {kTinnerf''^ that 
derives form the previous considerations reproduces quite 
well the overall increase of the gamma-ray luminosity of 
the ten clusters in Tab.l as a function of their inner core 
temperature kTinner- 



However, a scatter in the 



relation shown in 



FiglHl is indeed expected because the cluster gamma-ray 
luminosity depends on the different values of a and on the 
different values of the integration volume cx chosen for 
each specific cluster. The choice of different integration 
volumes for the considered clusters are motivated by the 
different extensions of the target protons in the thermal 
gas to produce effective — > tt" ^ 77 collisions, as in- 
dicated by the available X-ray data for each cluster (see 
Tab.l). A systematic increase of with Tinner is however 
shown by the data and can be understood by the previous 
scaling-law considerations. 

We have also calculated the gamma-ray flux expected 
from each cluster in Table [2] and we compared the results 
with the GLAST-LAT instrument sensitivity (see Fig.fTUl). 
While all the gamma-ray fluxes are lower than the avail- 
able EGRET hmits (Reimer et al. 2003), it is evident that 
6 clusters with Ky > 3 x 10~^ photons cm~^ s~^ (A262, 
A2199, Perseus, Hydra, A1795 and Coma), and 7 clusters 
with > 1.2 • 10~^ photons cm^^ s^^ (the previous clus- 
ters plus A2163) could be detected by the GLAST-LAT 
experiment in 1 year and 5 years at the 5a confidence 
level, respectively. 

The possible gamma-ray detection of these clusters will 
allow us to probe in details our model for the WR heating 
of the cluster cores and consequently set constraints on 
the parameters 7iwfl,o and a of the WR heating model. 

We stress that the gamma-ray emission evaluated here 
for cool-core clusters refer to the truly diffuse component 
arising from pp —+ —+ 77 decay. We do not consider here 
the possible contribution due to the gamma-ray emission 
from the central engine of the radio galaxies and / or from 
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Fig. 9. The correlation of the gamma-ray luminosity Lgamma 
with the inner temperature of the cool core Tinner - Blue dots 
refer to cool-core clusters and green dots to non cool-core clus- 
ters. 



their radio lobes. We will address more completely this 
issue in a forthcoming paper. 

Gamma-ray emission from WRs is able to probe di- 
rectly the hadronic component of their interaction with 
the IC gas. Nonetheless, the leptonic component of the 
WRs - IC gas interaction can be also efficiently probed 
by studying the secondary e^ e.m. emission: this is spread 
over a wide range of frequency, from radio (where the syn- 
chrotron emission dominates) to the UV, X-ray and HXRs 
(where the ICS emission dominates). We will discuss in the 
following these components separately. 

6. Diffuse radio emission in cool core clusters 

Diffuse radio emission is also inevitably produced by syn- 
chrotron emission of the secondary electrons produced by 
pp TT^ ^ interacting with the intra-cluster mag- 
netic field and it is therefore strongly related to the WR 
gamma-ray emission and to the temperature profile of the 
cluster core. From the WR spectrum and spatial distri- 
bution given in eq.®, we calculated the radio emission 
produced by the secondary electrons (we follow the pro- 
cedure described in Marchegiani et al. 2007) , by assuming 
a radial profile B{r) for the magnetic field. 

We apply this procedure to the clusters in Tab. [2] where 
a mini-halo is observed, i.e. A2390, Perseus, RX J1347.5- 
1145 and A2199 (in this last case an upper limit is only 
available, see Kempner & Sarazin 2000). 

In the cluster A2390, we found that the intensity of the 
radio emission is 2 x 10~^ and 2 mJy for values of central 
magnetic field in the range 1.3 - 13 (which are equal 
to i?eq and 10Bf,q, respectively, where Bf,q is the value in- 



Fig. 10. The distribution of the gamma-ray flux 
Fgamma{E > 0.1 GeV) with the inner temperature of the cool 
core -I inner • Blue dots refer to cool-core clusters and green dots 
to non cool-core clusters. We show the GLAST-LAT sensitivity 
for a 5(7 detection after 1 yr (cyan line) and after 5 yrs (green 
line) (from Ritz 2007, on behalf of the GLAST collaboration see 
http: / /www-glast .slac.stanford.edu / software/IS /glast _ lat-performance.htm 



ferred by equipartion between magnetic field and electrons 
energy, Govoni k Feretti 2004), and for a radial shape of 
the magnetic field B(r) cx ne{r), respectively. The radio 
flux associated to the A2390 mini-halo has been estimated 
to be 63 ± 3 mJy (Bacchi et al. 2003). However, the radio 
source at the center of A2390 is a complex combination of 
the radio emission from the powerful FR-II flat-spectrum 
cD radio galaxy and of the mini halo (Augusto, Edge and 
Chandler 2006). The substantial amount of polarization 
detected around the cD radio galaxy and in the mini halo 
region indicates that the radio flux attributed to the mini 
halo region could be substantially affected by the radio 
galaxy lobes, therefore pointing to a quite lower flux of 
the truly diffuse halo. 

For the Perseus cluster, we show in Fig. [11] the radio 
brightness at 1.4 GHz calculated for a magnetic field with 
central intensity of 10 /iC and a radial profile proportional 
to ng'^(r). These assumptions are comparable with those 
of Pfrommer & EnBlin (2004), which have performed a cal- 
culation of the radio emission of the secondary electrons in 
the Perseus cluster by fitting the CR density to obtain the 
radio intensity. We found that, for the same magnetic field 
they used [i.e. Bq = 10 ^G and B{r) oc n°-^{r)], the WR 
density needed to quench the CF to the observed temper- 
ature level reproduces very well the radio halo brightness 
(see Fig.lTTt. 

For the cluster RX J1347. 5-1145 we found that the 
shape of radio emission (see Gitti et al. 2007) can not be 
reproduced by secondary electrons if a monotonical radial 
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Fig. 11. Radio brightness at 1.4 GHz emitted by secondary 
electrons in Perseus cluster for Bq — 10 fiG and B{r) oc n2'^(r). 
The density and the radial shape of the WRs are those needed 
to stop the cooling flow at the observed temperature level. 



decrease of the magnetic field is assumed. Instead, if we 
assume Bq = I /xG and a profile of B{r) oc n~°-^(r), we 
obtain the brightness shape shown in Fig. [T^ This profile 
is very similar to that of Gitti et al. (2007), and the total 
flux calculated is 33 mJy, comparable to the measured one 
of ^ 25 mJy. It will be easy to recover the observed radio 
halo flux with a central magnetic field ^ 1 /iG. 

In A2199, where we use a radial profile B{r) oc nf.{r), 
we obtain a total fiux of 100 mJy for Bq = 5 /xG and 372 
mJy for Bq = 10 /xG. Since the upper limit is 168 mJy (see 
Kempner & Sarazin 2000), and these values of the central 
magnetic field are comparable with those obtained by FR 
measurements in cool-core clusters, we conclude that the 
available radio data are compatible with our model, if the 
central magnetic field is not very strong, i.e. for values of 
Bo ^7 fiG. 

For the two non cool-core clusters with radio ha- 
los considered in our analysis, we recover the results al- 
ready presented in previous analyses (Marchegiani, Perola 
& Colafrancesco 2007). Specifically, for Coma we found 
a — 1, i.e. the condition PwR oc Pth, and the same ratio 
Pwu/Pth ~ 0.3 (indicated also by the simulations of Ryu 
et al. 2003) for which the radio halo surface brightness 
can be recovered with a central magnetic field of 1.2 /xG 
(Marchegiani et al. 2007). 

For A2163 we again found a — I, i.e. the condition 
PwR oc Pth, and a ratio PwR/Pth ~ 0.2 for which the 
radio halo surface brightness can be recovered with a cen- 
tral magnetic field of 4.5 /xG, a factor ~ 2 higher than 
that found by Marchegiani et al. (2007). 



Fig. 12. Normalized radio brightness emitted by secondary 
electrons in the cluster RX J1347.5-1145 evaluated at 1.4 GHz 
for Bo = 1 /xG and B{r) oc n^'''^(r). The density and the radial 
shape of the WRs are those needed to stop the cooling flow at 
the observed level. 



We note that in clusters where no diffuse radio emis- 
sion has been detected, the magnetic field required by 
our WR model should be lower than the value required 
to make the radio emission from secondary electrons de- 
tectable. In fact, the synchrotron emission depends non- 
linearly on the magnetic field intensity and fluctuations 
(see discussion in Colafrancesco, Marchegiani & Perola 
2005) and therefore, little or moderate variations of its am- 
plitude can decrease the radio flux by even a large amount. 
It is important to remark that future, deeper radio obser- 
vation of galaxy clusters, as those which can be obtained 
with LOFAR and/or SKA, can provide further informa- 
tion on the intensity of the radio emission, and on the 
non-thermal content of galaxy cluster (see discussion in 
Marchegiani et al. 2007). 

We note also that the WRs model can provide a con- 
sistent picture for the formation of morphologically differ- 
ent radio sources in galaxy clusters: in fact, the diffusion 
processes in different clusters can produce different WRs 
radial proflles, then mini-halos or giant radio halos can 
be produced. In addition, the action of localized external 
compression (i.e. by action of merging) can enhance the 
CRs density and the magnetic field intensity in external 
regions of the clusters, by producing the radio relics. 
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7. Diffuse Hard X-Ray emission in cool core 
clusters 

The secondary electrons from WR interactions with the, 
IC gas protons also produce radiation in the Hard X-Ray 
band, by Inverse Compton Scattering (ICS) of the Cosmic 
Microwave Background (CMB) Radiation photons. 

The equilibrium spectrum of the secondary electrons 
depends on the density of the WRs (which is determined 
in our approach) and on the energy loss rates of the^ 
electrons (see, e.g., Marchegiani, Perola & Colafrancesco 
2007). While the ICS losses are known and are uniform all 
along the clusters volume, the synchrotron losses depend 
on i?^(r), and therefore they are not known, in princi- 
ple. However, for the clusters which host a diffuse radio 
emission, information on the magnetic field B(r) can be 
inferred from the observed radio-halo features (see Sect. 6 
above), while for the other clusters it is necessary to as- 
sume a magnetic field structure. 

Therefore, to predict the HXR emission of clusters 
listed in Tab. 1, we adopt two different strategies: i) for the 
radio emitting clusters we calculate the HXR emission by 
using a magnetic field compatible with the observed radio 
emission (see SectEl); ii) for all the clusters we calculate 
the HXR emission by assuming a low magnetic field as 
generally indicated by the ICS interpretation of the HXR 
emission observed in nearby clusters. In this last case, the 
synchrotron losses are everywhere weaker that the ICS 
ones, and the density of the secondary electrons - and 
consequently the HXR emission - have to be considered 
the maximum ones which is compatible with the cool-core 
heating by the WRs. The HXR fluxes reported in Tab. 2 
refer to this last evaluation strategy. 

In Figure [T3] the ICS emission spectrum of the sec- 
ondary electrons is shown for the radio emitting clusters 
(A2199, Coma, Perseus, A2163, A2390, RX J1347.5-1145). 
These fluxes - referring to the evaluation strategy i) - can 
be compared with the expected sensitivity in the HXR 
band of the next coming mission Simbol-X (Ferrando et al. 
2005), which is approximately ^ 10^^ photons cni^'^ s^^ 
keV~^ in the energy band 10—50 keV. From this Figure we 
can conclude that the clusters A2199, Coma and Perseus 
should be detectable by Simbol-X, if the magnetic fields 
in these clusters are those inferred by assuming that the 
radio emissions are produced by secondary electrons. The 
HXR emission expected in the 10-50 keV band should be 
clearly disentangled from the thermal one in this energy 
band. On the contrary, the clusters A2163, RX J1347.5- 
1145 and A2390 should not be detectable by Simbol-X in 
this evaluation strategy (for the case of A2390 we have 
used a magnetic field value ~ lOBeg, see discussion in 
Sect. 6). 

In Figlll] (both panels), we show the HXR ICS emis- 
sion of the secondary electrons in the WR model for all 
the clusters we considered (see Table [1]), by assuming a 
value of the magnetic field Bo = 0.1 //G, for which energy 
losses are completely dominated by ICS, i.e. our evalua- 
tion strategy ii). Note that these levels of HXR emission 
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Fig. 13. ICS emission in the HXR band of the secondary 
electrons for the radio emitting clusters. RX stands for RX 
J1347. 5-1145. The horizontal, long-dashed line represents the 
approximate sensivity of Simbol-X in the 10 — 50 keV range. 
The HXR fluxes have been computed under strategy i) (see 
text for details). 



should be considered as upper limits. From this figure, 
we can conclude that only RX J1347.5-1145 and A2390 
should be not detectable by Simbol-X; a difficult detec- 
tion is possible for A133 and A2163, while the other six 
clusters should be definitely detectable by Simbol-X, if 
the extreme conditions on the magnetic field that we have 
assumed hold. 

The HXR fiuxes predicted in the WR model are also 
much lower than the 20-80 keV fluxes of the nearby clus- 
ters detected with BeppoSAX (see Nevalainen et al. 2004). 
Specifically, we obtain for Coma a flux i^20-80fceV = 1-6 • 
2Q-13 gj.g g-i |-yg_ g^jj observed flux of F2o-80feey = 



'^■^ erg cm ^ s ^) 



(1.1 ± 0.5) • 10- 

F20-80keV — 2.1 

flux of F2o-80fcey = (8.4 ± 3.9) 
for A2163 a flux Fw-HOk.v = 3.2 • lO^^^^ 



10 erg cm ^ s ^ 



for A2199 a flux 
(vs. an observed 



10 erg cm ^ s ^) and 



erg cm ^ s ^ (vs. 
an observed flux of i^20-80fcev = (1.7±3.1)-10^^^ erg cm^^ 
s"^). The ICS HXR fluxes yielded by the WRs in cluster 
cores cannot reproduce the HXR fluxes observed towards 
these clusters with BeppoSAX (see also Colafrancesco, 
Marchegiani & Perola 2005 for a general, phenomenologi- 
cal discussion on this point). 

We also evaluated the non-thermal flux produced by 
ICS of the secondary electrons in the 2-10 keV energy band 
intergrated in a ~ 250 arcsec region around the center of 
Perseus to be compared with the thermal bremsstrahlung 
emission of the IC gas in the same region of the clus- 
ter. The non-thermal ICS flux induced by WRs is a factor 
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Fig. 14. ICS emission in the HXR band of the secondary elec- 
trons for all the clusters considered. For clarity, the emissions 
are reported in two different panels. RX stands for RX J1347.5- 
1145. The horizontal, long-dashed line represents the approx- 
imate sensitivity of Simbol-X in the 10 — 50 keV range. The 
HXR fluxes have been computed here under strategy ii) (see 
text for details). 



10^ less than the thermal bremsstrahlung flux. This fur- 
ther confirms that ICS HXR fluxes yielded by the WRs in 
cluster cores cannot reproduce the non-thermal flux ob- 
served towards the centers of clusters. 



We have shown that a specific feature of the WR scenario 
for the heating of cool cores in galaxy clusters is that it 
is inevitably associated to emission features that manifest 
themselves over a wide range of frequencies, from radio to 
gamma-rays. 

This fact provides a clear test to distinguish our model 
from other models of cooling-fiow quenching of mechani- 
cal and/or non- gravitational origin as well as from other 
models for the origin of non-thermal phenomena in the 
inner regions of galaxy clusters. 

How the WR model presented here compares with 
these models? 

Guo & Oh (2007) proposed a model of AGN feedback 
heating, in which cosmic rays produced by accretion trig- 
gered AGN activity heat the IC gas efficiently, with small 
dynamical perturbations on the IC gas itself. This model 
considers that the ICM is heated by cosmic-rays, which are 
injected into the ICM through diffusion or the shredding 
of the bubbles by Rayleigh- Taylor or Kelvin-Helmholtz 
instabilities, with the inclusion of thermal conduction. It 
seems that the inclusion of CRs in this model solves the ex- 
treme fine tuning required by the conduction-suppression 
factor when only thermal conduction operates, and pro- 
vides a more stable radial distribution of IC gas tempera- 
ture that is, therefore, mainly driven by the effect of CRs 
which set some low value of the mass accretion rate with 
little additional effect of the thermal conduction (see their 
Fig. 3 and their discussion). 

The large impact of CRs as the main heating and 
regulating agent of the clusters' cool-cores is quite sim- 
ilar to what is found in our model, because the con- 
duction effects are quite marginal. However, our WR 
model requires quite substantial WR pressures at the clus- 
ter's center (that becomes negligible in the outer regions) 
that are able to recover the non-thermal pressure pro- 
file found in Perseus and help reproducing the mini radio 
halo surface brightness of cool-core clusters (Perseus, RX 
J1347. 5-1145, A2199) with reasonable values and radial 
profiles of the cluster's magnetic field. The low CR pres- 
sure {PcR ~ O.lFt/i) and energy required by the model 
of Guo & Oh (2007) in A2199 can be explained as an 
effect of the lower energy cutoff {E„iin w 706 MeV) in 
the CR spectrum which is much higher than our choice 
{Emin = rupc^ ■ [1 + 3.4 • 10-^{kT/keV)] (see Sect.3). We 
have quantitatively shown this effect in the specific case 
of A2199 (see Fig. 6) where we compare our WR energy 
density distribution with that of the CR feedback heating 
model of Guo & Oh (2007). Such difference in the energy 
spectnim of CRs causes also a different level of predicted 
gamma-ray fiuxes for clusters with cool-cores, with our 
WR model providing higher fluxes detectable with the 
GLAST-LAT experiment and with the next generation 
HXR experiments like Simbol-X. 

In particular, the Perseus cluster should be one of the 
best candidates to look for a HXR ICS emission feature 
in the 10-50 keV range emerging from its inner core (see 
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Figs. 13 and 14). This prediction seems quite robust since 

the WR distribution responsible for such ICS emission is 
the same that is able to reproduce the thermal gas temper- 
ature distribution in its cool core and the mini radio halo 
surface brightness visible in the same spatial region. This 
prediction is also strengthened by the ability of the WR 
model to recover at the same time the thermal gas tem- 
perature radial distribution and the non-thermal (WR) 
pressure radial profile observed in the Perseus core (see 
Sanders & Fabian 2007). 

Other clusters with radio halos (e.g. Coma and the 
marginal detection of A2199) are also predicted to be vis- 
ible by Simbol-X once we normalize our model to their 
temperature profile and radio-halo emission. 

WR (CR) heating models are a viable alternative 
to mechanical heating models (e.g. Briiggen & Kaiser 
2002; Ruszkowski et al. 2004; Vernaleo & Reynolds 2006; 
Reynolds et al. 2005, McCarthy et al. 2007) in which the 
ICM is heated by the pdV mechanical work of the expand- 
ing bubbles, viscous dissipation of emitted soimd waves 
or mixing of the hot bubble plasma with the ICM. The 
detailed microphysics of how the mechanical heating pro- 
cesses take place has not been hammered out in detail, 
leaving a good deal of uncertainty; a definitive explanation 
for how energy is transported from the observed bubbles 
to the ICM in a distributed and isotropic fashion is still 
outstanding. The possible role of merging shocks to the 
heating of the IC gas in the clusters' cores is still to be 
quantitatively determined (see the simulations of Ryu et 
al. 2003). 

It is clear, nonetheless, that also cosmic-ray heating 
models suffer from other imcertainties: the details of how 
WRs adjust to a density profile similar to that of the IC 
gas have to be addressed by specific studies that could 
bring further robustness to the preliminary calculations 
presented here (see Appendix A and Sect. 2). 
It is possible to envisage various scenarios for the WR 
origin in terms of the activity of AGNs with small duty 
cycles and/or a substantial fractions of AGNs in cluster 
cores (see, e.g.. Bird et al. 2007), or CR origin in non- 
standard scenarios like cannonball models (sec deRi'ijula 
2004 for review), but we do not want to address this spe- 
cific topic here. However, it seems that the WR model 
is somewhat less uncertain that the CR feedback heat- 
ing, for which large uncertainties stand on how cosmic 
rays leak from the bubbles, and/or the rate at which bub- 
bles arc disrupted. It is nonetheless important to stress 
that elucidating the details of bubble disruption/cosmic- 
ray diffusion would be very useful in determining whether 
cosmic-rays or mechanical proccssscs provide a more effi- 
cient means of transporting heat from the bubble to the 
ICM. 

All of these issues may assume a great relevance in the 
light of the experimental possibility to study CR physics in 
large-scale structures in a multi-frequency approach: the 
next coming gamma-ray observations of several nearby 
clusters with GLAST, the future HXR detection of the 
ICS emission from the same clusters with Simbol-X and 



the possible detection of non-thermal SZ effect signatures 
(see, e.g. Colafrancesco, Marchegiani & Palladino 2003, 
Colafrancesco 2007 for a review) in the spectra of these 
clusters with SPT will clearly help to disentangle not only 
the presence of cosmic rays in the cluster atmospheres, 
but also the details of various WR models for the heating 
of cluster cores. 

In fact, models of mechanical heating do not show any 
non-thermal emission feature while the WR model pre- 
dict a substantial level of non-thermal emission at various 
frequency (radio, microwave, HXR and gamma) and the 
CR feedback heating model of Guo & Oh (2007) predicts 
only marginally visible non-thermal phenomena in cool- 
core clusters. 

The next generation radio, HXR and gamma-ray experi- 
ments will certainly disentangle these two classes of mod- 
els and will determine the amount of WRs which are 
present in galaxy clusters, provided that a good separation 
of the diffuse and point-like non-thermal emission will be 
achievable. 



9. Conclusions 

The WR heating model that we presented in this paper 
is able to reproduce the temperature, the pressure and 
the entropy structure of both cool-core and non cool-core 
clusters and it provides results that are consistent with all 
the available evidence of non-thermal phenomena emerg- 
ing from the cores of clusters, e.g. diffuse radio emission 
in the form of mini halos and diffuse radio halos, HXR 
emission limits obtained with BeppoSAX, Chandra and 
XMM, and gamma-ray limits obtained with EGRET. 

This model provides a theoretical description of the 
physics of cool cores that is directly related to the presence 
of several observable consequences of the presence of WRs 
in clusters' atmospheres. The main conclusions of our work 
are: 

— the presence of WRs in cluster atmosphere provides 
a simple solution for the temperature, pressure and 
entropy profiles of the thermal IC plasma and of the 
non-thermal plasma observed in the inner regions of 
several clusters studied in this paper. 

— The WR distribution produces a substantial emission 
of gamma-rays by their hadronic — > 7r° — > 77 inter- 
ations with the IC gas. A large fraction (6 clusters out 
of 10 considered in this study) have gamma-ray fluxes 
at £ > 100 MeV detectable with the GLAST-LAT ex- 
periment. The gamma-ray luminosity of the cool-core 
and non-cool core clusters correlates with the cluster's 
inner temperature, and thus provides a specific scaling 
behaviour of our WR model that can be verified with 
the next coming GLAST observations. 

— The WRs distribution can reproduce the radio emis- 
sion of all the cool-core clusters with mini-halos that 
we studied (A2199, Perseus, RX J1347.5-1145) except 
for A2390 (whose mini radio-halo flux could be affected 
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by residual radio emission from the lobes of the central 
cD radio galaxy). 

— The WRs distribution produces ICS emission with 
fluxes well below the HXR limits of A2199, Coma and 
A2163 provided by BeppoSAX and well below the not- 
thermal emission detected by Chandra and XMM in 
the central regions of Perseus. This means that ICS 
emission from WRs is not the explanation of the emis- 
sion of these clusters in excess with respect to their 
thermal bremsstrahlung radiation. 

— The distribution of WRs in clusters directly relates the 
thermal IC gas temperature, pressure and entropy dis- 
tribution in cool-core clusters (observable in the X-ray 
energy band) with non-thermal diffuse emission fea- 
tures observable at radio (synchrotron), HXRs (ICS), 
gamma-rays (mainly from tt'^ — > 77 decay) and mi- 
crowaves (by Sunyaev -Zel'dovich effect). We predict 
that the expected levels of these non-thermal emission 
features in the WR model will be testable with the 
next-coming experiments in the HXR band (Simbol- 
X), in the gamma-rays (GLAST-LAT) and in the mi- 
crowaves (SPT). These experiments can determine the 
amount of WRs which are present in galaxy clusters. 

— These observable predictions makes the WR model en- 
tirely testable by using a multi-frequency observational 
strategy that links, therefore, thermal and non-thermal 
phenomena in clusters cores. 

— The specific theoretical and observational features of 
the WR model render it quite different from other 
models so far proposed for the heating of cool-cores in 
galaxy clusters. Such peculiar difi^erences make it possi- 
ble to prove or disprove the WR model and, in general, 
the properties of any model that has been proposed as 
an explanation of the cooling-flow problems. 

Acknowledgements. We thank the Referee of this paper for var- 
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the manuscript. 
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Appendix A: CR proton diffusion in cluster 
atmospheres 

We calculate the effect of the diffusion of CR protons (we 
remind the reader that we use the same notation for CRs 
and WRs) in the cluster atmosphere according to the fol- 
lowing equation: 



where N(Ep,r) is the proton CR density, D{Ep,r) is the 
diffusion coefficient, bp(Ep,r) is the proton energy loss 
term and Qp{Ep, r) is the source term. 

Protons we are interested in here do not loose apprecia- 
bly their energy on time scales of the order of the cluster's 
lifetime so that we can neglect the term d{bpN)/dE in 
eq. (jA.l|l . We are interested in a quasi-stationary solution 
{dN/dt — 0), so that the diffusion equation simplifies to 



which, in spherical simmetry, can be written as: 
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The solution of this equation writes as: 
■ - / drr^Qp' 
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We search for a radial shape of the source term Qp{Ep, r) 
which, by the diffusion effect, gives a radial shape of the 
proton equilibrium distribution which is proportional to 
[g{r)]" (this is the radial dependence of the WRs to be 
inserted in eq.(|4])), where g{r) is the radial profile of the 
gas. Therefore, we write: 



r^D 
where 

/i(r) = 



drr'^Qp{r). 
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This is a general solution of eq. ljA.ip under the assumed 
hypothesis. Under the specific assumption that the diffu- 
sion coefficient is constant with the radius, we can write 

TT f dr'^^^NoWr. (A.7) 
-Do Jo ^ 

By deriving both members with respect to r, we obtain: 
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Do r 
and then 

h{r) = NoDoah{r), 
where 

h{r)^r'[g{rr-'^g{r). 

dr 

By deriving with respect to r, we obtain 
d 



Mr) 



-NoDoa—h{r) 
ar 



(A.; 



(A.9) 



(A.IO) 



(A.ll) 



and then we obtain the source term that satisfies eq. (jA.3|) : 



Qp{r) ^ -NoDoa\-^h{r). 



(A.12) 



(A.l) 



This equation gives the source term Qp{Ep,r) which pro- 
duces a radial distribution of protons proportional to 
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[(7(r)]", for a constant diffusion coefficient D{r) = Dq. In 
conclusion, we have shown that it is possible, by knowing 
the radial shape of the thermal gas in a cluster, to know 
the radial shape that the WR source term should have 
in order to recover the required radial distribution of the 
warming rays. 

In the following, we consider a solution of eq. (jA.3P 
for a specific case of g{r), and we find the correspond- 
ing shape of the source term Qp{r). We consider only the 
radial shape of the proton distribution, and we use a free 
normalization. 

We consider a thermal IC gas profile of the form 



air) 



(A.13) 



Equations (|A.12p and (jA.lOp provide the shape of the 
source term Qp{r). In Fig lA.H we consider the case of a 
cluster with rc = 0.3 Mpc and qth = 1.125 (corresponding 
to the case of Coma). As we can see from this figure, the 
effect of the proton diffusion is to provide a broad equi- 
librium distribution of CRs from a narrow distribution of 
the proton source term. 



with a small value of is able to reproduce the source 
term required to provide the CR profile [^(r)]" with a = 1. 
In Fig lA.2| we compare the source term in Fig lA.ll with a 
Gaussian distribution as in ea. (jA.14p with — 0.17 Mpc 
and two other Gaussians with ~ 0.10 and = 0.05 
Mpc. Then this function, with properly chosen parameters 
(normalization and rg), can be a good approximation for 
the required source term. 

In Fig. lA.3l we show the equilibrium distribution of the 
GR protons produced by several Gaussian source terms as 
in eq. (|A.14|) with = 0.17, 0.10 and 0.05 Mpc, and we 
find that the first function is very similar to the IC gas ra- 
dial distribution of a Coma-like cluster. Since this proce- 
dure is not an accurate fit, but only a qualitative compar- 
ison for an arbitrary set of parameters, we can conclude 
that a narrow CR proton source term, by effect of the 
spatial diffusion, can produce a broad radial equilibrium 
distribution of the CR protons which is very close to the 
thermal IC gas distribution. 
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Fig. A. 2. Source term (dot-dashed line) as in Fig lA.ll com 

pared with a Gaussian shape as in ea. (|A.14[) with 
Mpc (solid line), — 0.10 Mpc (dashed line), and r, 
Mpc (dotted line). 



0.17 
0.05 



Fig. A.l. Source term [ea. HA.12|) ] for the case of a cluster 
with the thermal profile as Coma {rc — 0.3 Mpc and qth = 
1.125) and for a = 1. With the dashed line, it is drawn the 
thermal profile of the cluster (which, in this case, is also the 
resulting proton radial distribution), with an arbitrally chosen 
normalization for a better comparison. 



A Gaussian distribution 

2' 



j{r) (X exp 



(A.14) 
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10. OOF 




Fig. A. 3. Equilibrium distribution of the protons produced 
by a Gaussian source term as in eq. HA.14|) with = 0.17 Mpc 
(solid line), Vs — 0.10 Mpc (dashed line), and rs = 0.05 Mpc 
(dotted line). 



